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STABILITY ANALYSIS 
SLOPES AND EMBANKMENT FOUNDATIONS 


PROGRAM NUMBER 5171 


ABSTRACT 


THE COMPUTER PROGRAM PRESENTED HERE CAN BE USED WITH THE 1BM 
1620 COMPUTER TO ANALYZE THE STABILITY OF CUT SLOPES, FILL SLOPES, AND 
EMBANKMENT FOUNDATIONS. CROSS=8ECTIONS CONTAINING ANY IRREGULAR CON@= 
FIGURATION OF GROUND SURFACE AND SUBSO}L LAYERS THAT CAN 8E APPROXIMATEO 
BY TEN STRAIGHT LINES CAN BE ANALYZEO. A FAILURE ARC ANDO A FACTOR OF 
SAFETY SUFFICIENTLY CLOSE FOR PRACTICAL PURPOSES TO THE MOST CRITICAL 
FAILURE ARC ANO THE LOWEST FACTOR OF SAFETY ARE FOUND AUTOMATICALLY. THE 
COMPUTER TIME REQUIREO FOR A TYPICAL PROBLEM 18 40 MINUTES.1F THE APPROXI= 
MATE LOCATION OF THE MOST CRITICAL ARC CAN BE PREDICTED. OTHERWISE, THE 
TIME REQUIREMENT MAY BE CONSIDERABLY GREATER. 

AS WITH ANY STABILITY ANALYSI8, THE UNDERSTANDING OF THE UNDER= 
LYING ASSUMPTIONS AND THE CHOICE OF THE APPROPRIATE SHEAR STRENGTH PARA= 
METERS ARE ESSENTIAL TO THE OBTAINING OF A REALISTIC ANSWER TO THE PROPER 
EVALUATION OF THE RESULTS. NO ATTEMPT HAS BEEN MADE IN THIS REPORT TO 


GiVE THE READER THE NECESSARY BACKGROUND KNOWLEOGE OF SOJ].L MECHANICS. 
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INTRODUCTION 
A NUMBER OF COMPUTER PROGRAMS HAVE BEEN DEVELOPED TO REDUCE THE 
TIME REQUIRED FOR THE STABILITY ANALYS18 OF SLOPES AND EMBANKMENT FOUNDA= 
TIONS. THE PROGRAM PRESENTED IN THIS REPORT 18 BASEO ON THE PROGRAM 
"Computer SovtuTIons FOR Store Stasicity ANALYSIS = 9.2.020" WRITTEN BY 
tHE Corps oF ENGINEERS, U.S. ARMY ENGINEER Dig TRICT, ALBUQUERQUE, NEw 
Mexico, FoR THE IBM 1620 compuTER. THE FOLLOWING TWO OPERATIONS, NOT 
POSSIBLE WITH THE ORIGINAL PROGRAM, HAVE BEEN INCLUDED IN THE PRESENT PROGRAM: 
1. THE COMPUTATION OF THE STABILITY OF AN EMBANKMENT ON A 
FOUNDATION CONTAINING ONE OR MORE INCOMPLETELY CONSOLI= 
DATED SOIL LAVERS (U 18 Less THAN 100%). 
2. AN AUTOMATIG SEARCH FOR THE MINIMUM FACTOR OF SAFETY. 
To 0O THIS WITHOUT OVERRUNNING THE AVAILABLE STORAGE CAPACITY AND 


TO SHORTEN THE COMPUTING TIME, THE ORIGINAL PROGRAM HAD TO BE MODIFIED 


CONSIDERABLY. 


Basis FOR PROGRAM 

THE SLOPE STABILITY ANALYS!8 USEO IN THE PROGRAM 18 THE CIRCULAR 
ARG METHOD USING FINITE SLICES AND CONTAINING THE SIMPLIFYING ASSUMPTION 
THAT THE RESULTANT OF FORGES ON THE SIDE OF A SLICE 18 ZERO. THIS ANALYSIB, 
ALTHOUGH CONSERVATIVE IN CASES WHERE THE ANGLE OF INTERNAL FRICTION IS 
GREATER THAN ZERO, HAS THE ADVANTAGE OF BEJNG MUCH LESS COMPLICATED THAN 
MORE RIGOROUS ANALYSES WHEN A NUMBER OF SOIL LAYERS WITH DIFFERENT STRENGTH 
CHARACTERISTICS ARE PRESENT. 

New York State DEPARTMENT OF PUBLIC WORKS USES THE TOTAL STRESS 
METHOO OF STABILITY ANALYSIS. THE ANGLE OF INTERNAL FRICTION AND THE CO= 
HESION FOR USE IN THE COMPUTATION OF STABILITY OF EMBANKMENTS ARE DETERMINED 
FROM THE TOTAL STRESS FAILURE ENVELOPE FOR CONSOLIDATED@“UNDRAINED TRIAXIAL 


TESTS. TO TAKE INTO ACCOUNT THE GRADUAL STRENGTH GAIN OF A SATURATED COHESIVE 
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FOUNDATION SO}L DUE TO CONSOLIDATION UNDER A SUPERIMPOSED LOAD, THE COM= 
PONENT OF THE LOAD NORMAL TO THE. FAILURE ARGC IN TH& FOUNDATION SOJL 18 
MULTIPLIED BY THE DEGREE OF CONSOLIDATION (U) OF THE 80/tL. THE 8UM OF THE 
NORMAL COMPONENT OUE TO THE SUPERIMPOSED LOAD TIMES THE DEGREE OF CONSOLI-= 
DATION ANDO THE NORMAL COMPONENT OUE TO PREEXISTING OVERBURDEN 18 MULTIPLIED 
BY THE ANGLE OF INTERNAL FRICTION TO GIVE THE FRICTIONAL RESISTANCE. 

THE DEGREE OF CONSOLIDATION IN SOIL MECHANICS 18 GENERALLY COMPUTED 
WITH RESPECT TO THE CHANGE IN VOID RATIOc IT MUST BE KEPT IN MIND, HOWEVER, 
THAT THE RELATIONSHIP BETWEEN VOID RATIO AND EFFECTIVE STRESS AND STRENGTH 
18 LOGARITHMIC. AT INTERMEDIATE DEGREES OF CONSOLIDATION BETWEEN O ano 1.0, 
THE DEGREE OF CONSOLIDATION WITH RESPECT TO THE CHANGE §N VOJDO RATIO CAN 
BE CONSIDERABLY HIGHER THAN THE DEGREE OF CONSOLIDATION WITH RESPECT TO THE 
CHANGE IN EFFECTIVE 8TRES8 AND STRENGTH, THEREFORE, THE TERM "DEGREE OF 
CONSOL! DATION (U)" as USED IN THI8 REPORT DESIGNATES THE RATIO OF THE EFFEC= 
TIVE STRESS INCREASE AT ANY TIME TO THE EFFECTIVE STRESS INCREASE AT FULL 
CONSOLIDATION. 

THE SEARCH PATTERN FOR THE MOST CRITICAL FAILURE ARC 18 BASED ON THE 
EXISTENCE OF "CONTOUR LINES" OF FACTORS OF SAFETY. POINTS ARE INVESTIGATED, 
GOING IN THE DIRECTION OF DECREASING FACTOR OF SAFETY, UNTIL A POINT 18 FOUND 
WHERE THE FACTOR OF SAFETY 18 LOWER THAN AT SURROUNDING “pig THE SEARCH 


PATTERN 18 DESCRIEED MORE COMPLETELY UNDER SEQUENCE OF COMPUTATIONS. 
CAPABILITIES AND LIMITATIONS OF PROGRAM 

THE PROGRAM CAN ANALYZE A CROSS=SECTICN CONTAINING INCLINED AND 
DISCONTINUOUS FOUNDATION LAYERS, AN IRREGULAR GROUND SURFACE, ANO BERMS. 
THE LIMITATIONS OF THE PROGRAM ARE AS FOLLOWS: 


ir 


ORN, J.A. “Computer ANALYSIS OF SLOPE STaBiILi TY" 
JOURNAL, SOILS MECHANICS AND FOUNDATIONS DIVISION, 
ASCE, Vor. 86, No. SM3, June 1960. 
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FoR USE wiTH THE Basic !BM 1620 computer (20,000 storace 
POSITIONS) THE NUMBER OF LINES MAKING UP THE CROS8=SECTION 
18 LIMITED TO TEN, ANO THE NUMBER OF INTERSECTIONS BETWEEN 
THE FAILURE ARC AND LINES MAKING UP THE CROSS<SECTICN 18 
LIMITED TO FOURTEEN. THI& MEANS THAT A CROSS=SECTION CON@= 
TAINING AN EMBANKMENT AND CONTINUOUS FOUNDATION LAYERS 
WOULD BE LIMITED TO SIX FOUNDATION LAYERS. IF A BERM WERE 
INCLUOED, THE NUMBER OF FOUNDATICN LAYERS PERMISSIELE WOULD 
BE REDUCED TO FIVE. 


No VERTICAL LINES OR LINES HAVING A SLOPE GREATER THAN 500 

ON | SHOULD BE USED. IF A VERTICAL LINE CONSTITUTES A PART 
OF THE FREE SURFACE OF THE SECTION, IT SHOULD BE APPROXIMATED 
BY A LINE WITH A STEEP SLOPE. OTHERWISE, A VERTICAL LINE 
BETWEEN DIFFERENT ZONES IN THE SECTION CAN SE LEFT OUT. 


THE MAXPMUM VALUES WHIGH MAY BE ASSIGNED TO X=-AND Y=COOROINATES 
ARE PLUS OR MINUS 9399. 


THE HORIZONTAL DISTANCE BETWEEN THE CENTER OF FAILURE ARC 
GIVEN IN THE INPUT ANDO THE ACTUAL MOST CRITICAL CENTER AS 
DETERMINED BY THE COMPUTER CANNOT BE GREATER THAN Forty (40) 
FEET. 


CURVED LINES MUST BE APPROXIMATED BY STRAIGHT LIMES. 


SEEPAGE FORCES ARE NOT TAKEN INTO ACCOUNT. HOWEVER, THEIR 
EFFECT CAN BE APPROXIMATEO BY USING IN THE SEEPAGE ZONE 
THE SATURATED UNIT WEIGHT AND A DEGREE OF CONSOLIDATED (U) 
EQUAL TO THE BUOYANT UNIT WEIGHT DIVIDEO BY THE SATURATED 
UNIT WEIGHT. AS EXPLAINED LATER, THE UNIT WEIGHT IN THIS 
CASE MUST BE WRITTEN ON THE INPUT FORM WITH A MINUS SIGN 
PRECEDING IT. 


FOR A HEIGHT OF SLOPE GREATER THAN TEN FEET THE DIFFERENCE BETWEEN 


THE MINIMUM FACTOR OF SAFETY FOUND BY INVESTIGATING POINTS FIVE FEET APART, 


AS IN THI& PROGRAM, AND THE MINIMUM FACTOR OF SAFETY FOUND BY INVESTIGATING 


POINTS ONE FOOT APART HAS BEEN FOUND DURING THE OEVELOPMENT OF THIS PROGRAM 


TO BE LESS THAN FIVE PERCENT. 


INPUT INSTRUCTIONS 


THE FOLLOWING STEPS SHOULD BE FOLLOWED WHEN PREPARING AN INPUT FORM 


FOR COMPUTER ANALYS183 


1. 


DRAW THE SECTION ON GROSS=SECTION PAPER 80 THAT THE CENTER 
OF THE MOST PROBABLE FAILURE ARC 18 TO THE RIGHT OF THE 
HIGHEST POINT ON THE GROUND SURFACE. 
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NUMBER ALL LINES MAKING UP THE SECTION CONSECUTIVELY, 
STARTING WITH THE TOPMOST LINE AS NUMGER ONE. No vUINe 
SHALL, AT ANY POINT, BE DIRECTLY ABOVE A LOWER=NUMBERED 
LINE. 


USING THE HIGHEST POINT ON THE GROUND SURFACE A& ORIGIN, 
ASSIGN COORDINATES TO THE ENOPOINTS OF, AND INTERSECTIONS 


BETWEEN, ALL LINES. X 18 POSITIVE IN THE OIRECTION FROM THE 


ORIGIN TOWARDS THE ASSUMED CEMTER OF THE FAILURE ARC, OR 
TO THE RIGHT. Y 18 POSITIVE UPWARDS. THE YeCOORDINATES 
OF ALL LINES MAKING UP THE SECTION ARE ZERO OR NEGATIVE. 


CHOOSE A STARTING FAILURE ARC. IN ORDER TO SAVE TIME IT 
1S OESIRABLE THAT THIS FAILURE ARC BE AS NEAR THE MOST 
CRITICAL ARC AS POSSIELE. HOWEVER, SINCE IN THE SEARCH 
PATTERN THE YeCOORDINATE OF THE CENTER OF FAILURE ARC CAN= 
NOT BE DECREASEG, THE STARTING CENTER SHOULD BE CHOSEN 
NOT HIGHER THAN THE LOWEST FPOSSIGLE POSITION OF THE MOST 
CRITICAL CENTER. ALSO, SINCE NO FAILURE ARC EXANINED 
DURIAG THE SEARCH PATTERN 18 SHALLOWER THAN THE STARTING 
FAILURE ARC, THE STARTING RADIUS SHOULO BE CHOSEN 60 
THAT THE STARTING FAILURE ARC 18 NOT BELOW THE PROBABLE 
POSITION OF THE MOST CRITICAL ARC. 


IN A SOIL DERIVING MOST OF ITS STRENGTH FROM INTERNAL 
FRICTION, A SHALLOW FAILURE ARC WILL BE THE MOST CRITICAL; 
IN A SOUL DERIVING MOST OF ITS STRENGTH FROM COHESION, A 


‘DEEP FAILURE ARC WILL USUALLY BE THE MOST CRITICAL. For 


THIE REASON, IN A CASE WHERE A COHESIVE MATERIAL 18 OVER] 
LAI" BY A COHESIONLESS MATERIAL, IF THE STARTING FAILURE 
ARC LIES IN THE COHESIONLE&SS LAYER, THE COMPUTER USING 
THIS SEARCH PATTERN MAY FINO THE MOST CRITICAL CENTER FOR 
FAILURE §N THAT LAYER WITHOUT INVESTIGATING THE COHESIVE 
LAYER. TO PREVENT THI8, JUOGMENT MUST BE EXERCISED IN THE 
SELECTION OF THE STARTING FAILURE ARC. 


IN SOME CASES, SUCH AS THAT OF AN EMBANKMENT WITH A 
STABILIZING BERM, MORE -THAN ONE “LOW SPOT" MAY EXIET ON 
THE SURFACE FORMED BY THE FACTOR OF SAFETY CONTOUR LINES, 
AND IT MAY BE DESIRABLE TO ANALYZE THE SECTION USING MORE 
THAN ONE INITIAL FAILURE ARC CENTER. 


IF A THIN, LOW“STRENGTH SOIL LAYER 18 PRESENT, THE 
STARTING RADIUS SHOULD BE SHOSEN SUCH THAT THE RADIUS PLUS 
A MULTIPLE OF FIVE FEET WILL FALL WITHIN THE LAYER. 


Fill OUT THE INPUT DATA SHEET. WRITE DOWN THE NUMBER OF LINES 


IN THE SECTION. List THE COORDINATES OF EACH LINE AND THE 
SO!L PROPERTIES IN THE SOIL BELOW EACH LINE IN THEIR RE= 


SPECTIVE COLUMNS AFTER THE CORRESPONDING LINE NUMBER, WRITE 
IN THE APPROPRIATE SPACES THE COORDINATES ANDO RADIUS DEFINING 


THE STARTING FAILURE ARC. FOLLOWING 18 AN EXPLANATION OF 
THE TERMS USED IN THE INPUT DATA SHEET? 
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NUMBER OF LINES 


X} = LEFT HORIZONTAL COORDINATE OF A LINE. 

Yj = LEFT VERTICAL COORDINATE OF A LINE. 

X2 = RIGHT HORIZONTAL COORDINATE OF A LINE 

Y2 = RIGHT VERTICAL COORDINATE OF A LINE. 

W = UNIT WEIGHT OF THE MATERIAL BELOW THE LINE. 
FOR 8OILO BELOW THE GROUND WATER TABLE, THE 
BUOYANT UNIT WE!GHT SHOULD BE USED. IF A 
FOUNDATICN LAYER 18 NOT FULLY CONSOLIDATED 
UNDER THE WEIGHT OF AN OVERLYING SO!fL LAYER, 
THE UNIT WEIGHT OF THE OVERLYING SOIL LAYER 
18 WRITTEN WITH A MINUS SIGN IN FRONT OF IT. 

F = TANGENT OF THE ANGLE OF INTERNAL FRICTION 
OF THE MATERIAL BELOW THE LINE.) 

C = COHESION OF THE MATERIAL BELOW THE LINE. 

U = OEGREE OF CONSOLIDATICN, EXPRESSED AS A 
DECIMAL, IN THE MATERIAL BELOW THE LINE. 

THIS VALUE MUST BE FILLED JN FOR ALL 80IL8 
IN THE SECTION. 

H = HORIZONTAL COORDINATE OF THE STARTING FAILURE 
ARG CENTER. 

G = VERTICAL COORDINATE OF THE STARTING FAILURE 
ARC CENTER. 

R = STARTING RADIUS. 


Ve NEGATIVE NUMBERS MUST BE WRITTEN WITH A MINUS SIGN PRECEDING 
THEMS POSITIVE NUMBERS SHOULD BE WRITTEN WITH NO SIGN. 


8. ALL INPUT DATA, EXCEPT THE NUMBER OF LINES (N) must BE 
WRITTEN WITH A DECIMAL POINT, AS SHOWN ON THE SAMPLE INPUT 
SHEET. 


SEQUENCE oF CompuTATIONS 


THE GENERAL SEQUENCE OF STEPS FOLLOWED BY THE COMPUTER IN FINOING 
THE MOST CRITICAL FAILURE ARC AND THE MINIMUM FACTOR OF SAFETY FOR A SLOPE 
ARE AS FOLLOWS: 


1. THE INTERSECTION OF EACH LINE DEFINED BY THE ENDPOINT 
X ANO Y COORDINATES RECEIVED ON THE INPUT CARDS, WITH 


THE FAILURE ARC, DEFINED BY THE HORIZONTAL AND VERTICAL 
COORDINATES OF THE CENTER ANO BY THE RADIUS (H, G ano R), 


§8 FOUND. 
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THE INTERSECTIONS ARE ARRANGED IN ORDER OF DECREASING 
MAGNITUDE OF X=COORDINATE ANO NUMBERED. 


EACH SEGMENT OF ARC WIDER THAN ONE=TENTH OF THE RAOIUS 
8&8 DIVIDED INTO TEN “SLICES” OF EQUAL WIDTH. EACH 
SEGMENT OF ARC EQUAL TO OR SMALLER IN WIDTH THAN ONE= 
TENTH OF THE RADIUS 18 DIVIDEO INTO TWO “"sL;3CE8" OF 


“EQUAL WIDTHe 


THE Y=COORDINATES OF THE INTERSECTIONS OF THE CENTER= 
LINE OF THE SLICE ANDO THE LIWES MAKING UP THE SOIL SEC= 
TION ARE DETERMINED.) 


THE WEIGHT OF EACH SLICE TO BE USED JN COMPUTING THE 
OVERTURNING FORCE 18 COMPUTED BY MULTIPLYING THE Y= 
DISTANCE BETWEEN SUCCESSIVE INTERSECTIONS ON THE CENTERLINE 
BY THE SLICE WIDTH AND BY THE UNIT WEIGHT OF THE CORRES= 
PONDING MATERIAL ANO ADDING THE RESULTS FOR ALL PORTIONS 

OF THE CENTERLINE. 


THE WEIGHT OF EACH SLICE TO BE USED IN COMPUTING THE 

NORMAL FORCE ON THE ARC 18 COMPUTED. IN THE CASE OF 

COMPLETE CONSOLIDATION THIS WILL BE EQUAL TO THE WEIGHT 
COMPUTED FOR OVERTURNING FORCE. IF THE FAJLURE ARE PASSES 
THROUGH A LAYER THAT 18 NOT FULLY CONSOLIDATED, THE WEIGHT 

OF THE PORTION OF THE SLICE PASSING THROUGH THE LAYER OF 80IL 
UNOER THE ACTION OF WHICH INCOMPLETE CONSOLIDATION EXISTS 

18 MULTIPLIED BY THE DEGREE OF CONSOLIDATION (U). 


THE COHESIVE FORCE ACTING ON THE BASE OF A SLICE 18 COM= 
PUTED BY MULTIPLYING THE UNIT COHESION (C) BY THE CHORD 
LENGTH OF THE BASE OF THE SLICE. THE ERROR INTROOUCED BY 
USING THE CHORD LENGTH INSTEAD OF THE ARC LENGTH 1S LESS 
THAN ONE PERCENT. 


THE OVERTURNING AND THE NORMAL COMPONENTS OF THE WEIGHT OF 
EACH SLICE ARE OETERMINEO. THE FRICTION FORCE 18 COMPUTED 
BY MULTIPLYING THE NORMAL COMPONENT BY THE ANGLE OF INTER@= 
NAL FRICTION, 


AS THE COMPUTATIONS PROGRESS FROM SLICE TO SLICE, THE TOTAL 
OVERTURNING, FRICTION, AND THE COHESIVE FORCES ARE FOUND BY 
ADOING THEM CUMULATIVELY FROM 8LICE TO SLICE, 


THE FACTOR OF SAFETY 18 COMPUTED BY DIVIDING THE SUM OF 
THE FRICTION ANDO THE COHES!VE FORCES BY THE OVERTURNING 
FORCE. 


THE RADIUS 18 INCREASED BY FIVE FEET AND STEPS ONE THROUGH 
ELEVEN ARE REPEATED TO FIND THE FACTOR OF SAFETY FOR THE 
NEW RADIUS. THiS 18 REPEATED AS LONG AS EACH SUCCEEDING 
FACTOR OF SAFETY 18 LOWER THAN THE PREVIOUS ONE. WHEN A 
FACTOR OF SAFETY LARGER THAN THE PRECEDING ONE 1S FOUND » 
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THE PRECEOING FACTOR OF SAFETY ANO THE X= ANDO Yeo 
COORDINATES AND THE RADIUS OF THE FAILURE ARC ARE PUNCHED 
ouTe NEXT, THE X=COOROINATE OF THE FAILURE ARC CENTER 18 
INCREASED BY FIVE FEET, THE RADIUS 18 SET EQUAL TO THE INPUT 
RADIUS ANO THE LOWEST SAFETH FACTOR FOR THESE CENTER COOR= 
DINATES DETERMINED ANDO PUNCHED OUT. ADDITIONAL FAILURE ARC 
CENTERS ON A HORJZONTAL LINE ARE EXAMINEO UNTIL THE MOST 
CRITICAL CENTER ON THAT LINE 18 FOUND. THEN, ITS X= 
COORDINATE 15 HELD CONSTANT, BIT ITS YeCOORDINATE }NCREASEO 
BY FIVE FEET AND THE MOST CRITICAL CENTER ON THAT HORI ZON= 
TAL LINE 18 FOUND. THIS PROCESS 18 CONTINUED UNTIL A 

POINT 1§8 FOUND AT WHICH, IF IT 18 ABOVE THE INPUT POINT, 
THE FACTOR OF SAFETY 18 LOWER THAN AT ALL POINTS SURROUND= 
1 G §T ON A FIVE FOOT GRID, HORIZONTALLY, VERTICALLY AND 
DIAGONALLY, EXCEPT THE GRID POINT FIVE FEET ABOVE AND FIVE 
FEET YO THE LEFT OF ITe 
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SLOPE STABILITY ANALYSIS 


PROGRAM NUMBER #5171 


SAMPLE PROBLEM 


Jos TRANSMITTAL Caro 

Cross Section DiAGRam 

Source SHeet - SM=-285 (9/63) 
Output LISTING 


Output Caros 


OCTOBER 1963 
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SAMPLE JOB TRANSMITTAL FORM 


SLOPE STABILITY ANALYSIS 


PROGRAM # 5171 








Form EDP - 3a (7/63) 
BUREAU OF EDP 


JOB TRANSMITTAL CARD | Control Identification Number 


Submitting Bureau or Unit SOs MECHAN! g 5 Location i PAD So ee 
Oslo GiB NM wivoe § topes 2... 2404. 7 


Date of Submission 
Date of Completion _ol9loielaial | Actual Completion Date | | | | | | | 


Job Title: , SLOPE ef LAB Lf / Y ANAL 1S Vs) Procedure No. 


Remarks: 
Recurring 
[] Special. 
TX key Punch p= Tabulate - CJ Reproduce Ea Sort 
C] Printing Punch C] Interpret ie Collate AL compute 


For further detailed instructions, see reverse side. 


SEE REVERSE SIDE OF THIS FORM FOR DETAILED INSTRUCTIONS 
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STATE OF NEW YORK DEPARTMENT OF PUBLIC WORKS BUREAU OF SOIL MECHANICS 


PROJECT___ SAMPLE PROBLEM 


STATION 
ASSUMED CONDITIONS 


SLOPE STA 
PROGRAM NO. SI7I 


CARD TYPE 2 











Prepared by R. GODELL pate-O/3/63 


BILITY ANALYSIS 
1.B.M. |620 COMPUTER 
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OUTPUT CARDS FOR SLOPE STABILITY ANALYSIS PROGRAM 5171 





y, fy 2 3 4[5 6 7 8[9 10 11 12|13 14 15 16M 16 19 20[21 22 29 24[25 25 mM 20/29 30 31 32]33 34 35 pR]37 3a 39 40]47 42 40 44] 46 47 40[49 50 51 52]59 54 55 56[57 50 59 G0]G1 62 69 64]65 66 G7 68/69 70 71 72[73 74 75 76]77 78 79 80] 


a a 
| 000 0/0 00/000 0/0000/00 0000 0/0008 000/000 0MM 00/000 00MM 0/000 0]0000/0000/0000\0000/000000008Hm0 
1234/5 67 
Ti1igit 


819 10 11 12]13 14 15 16/17 18 19 20/21 22 23 24]25 26 27 26/29 30 31 32/33 34 35 36/37 38 39 40/41 42 43 44]45 46 47 48]49 50 51 52/53 54 55 56/57 58 59 60/61 62 63 64/65 66 67 68/69 70 71 72173 74 75 76)77 78 79 80 


| Tt PRO eRe PRP POR MOR CRO ORO | PPR CR Oe Pee Ree ee eee eee etme | 
| 2222/222 21222 2122 2 2/222 2/22 2 212 2 2 212 2 2 22 22 2/2 2 2 2\2 2 2 2/2 2 2 22 2 2 242 2 2 212 2 2 212 2-2 22 2 2 22 22 2/2 2 2 2M 2 2 2 


12.3.415 6 7 819 10 19 12}13 14 15 16/17 18 19 20/21 22 23 24/25 26 27 28/29 30 31 32/33 34 35 36}37 38 39 40)41 42 43 44145 46 47 48/49 50 51 52/53 54 55 56/57 58 59 60/61 62 63 64/65 66 67 68/69 70 71 72)73 74 75 76|77 78 79 80 


330) 3/333 3/3 33 3/33.339 333/333 39 3M 3/3 33 3/3 3 38/3 3 3 3/3 33 3h 3 3 3/3 3 3 3)3 3 3 3/3 3 3 3/3 3 3 3/33 3 3/3 3 3 3/3 3 3:33.33 3 
4444)444414444)444 4184441444 414444/444 4194441444 4)444044 4 41444 4144441444 41444 414 4 44 4 4 4444 4 4d 


12 3 415 6 7 BIS 10 19 12/13 14 15 16/17 18 19 20)21 22 23 24125 26 27 28/29 30 31 32)33 34 35 36/37 38 39 40/41 42 43 44145 46 47 48/49 SO 51 52/53 54 55 56/57 58 59 60/61 62 63 64/65 66 G7 68/69 70 71 72|73 74 75 76)77 78 79 80 


§ 555/555 51555 515 5 5 515 5 6 5/5 5 5 5/5 5 5 5/5 5 5 515 6 5 515 5 6 515 6 5 515 5 5 5/5 5 5 5/5 5 5 5/55 5 9/555 9/55 5 5/55 5 9/555 55555 
GHCEG GEG GE GGG EGG G6 GGECCGG GGG E66 CECE CGE GEC CCE C6666 ESIC C CCG 6666666666 6666 666666666 


12.39 4/5 6 7 BI9 10 11 12]13 14 15 16}17 18 19 20)21 22 23 24)25 26 27 26/29 30 31 32/33 34 35 36)37 38 39 40/41 42 43 44/45 46 47 48/49 50 51 52/53 54 55 56157 58 59 60/61 62 63 64/65 66 67 63/69 70 71 72/73 74 75 76\77 78 79 BO 


| 17790779/7779)77790)777.7)777: 17 7 777000 0) 700) 7 0 00 77 77 7770777 





818.88 8|8 88 CH MH 8/8 8 8 8/8 GM bis 6 8 8/5 8 OM i3 8 8 8/8 88 eM SB 8i8 f 8 Bik GS G/B 8 Bl8 8 8 8 8 8 8/8 BB Bi 8 BC 8 8B 


49 50 51 52/53 54 55 56/57 58 59 68/61 62 63 64/65 66 67 68/69 70 71 72/73 74 75 76)77 78 79 80 


99991999 $9959 9199991999 91999 91999919999 


49 50 51 52153 54 55 56/97 58 59 66161 62 63 64165 66 67 68169 70 71 72173 74 75 16177 78 79 85 


17 18 19 20/21 22 23 24)25 25 27 28/29 30 31 32133 34 35 36/37 38 39 46 


999919999199 99/9 9:9 91999919999 


45 46 47 46 


Nee! 


45 46 47 43 


41 42 43 44 


993.9 


41 42 43 44 


8/9 10 11 12)13 14 15 16 


8 8 

67 

999/999 919999 
8 


























GENERAL PURPOSE, 20 FIELD 


1 B19 10 11 12]13 14 15 16117 18 19 20121 22 23 24125 26 27 28129 30 31 32133 34 35 35137 38 39 40 
| BM 733727 ; 


{12 3 4]5 6 7 8]9 fo 11 12]19 14 15 16 18 19 20] 18 19 20]21 22 23 24]25 25 BY 26]29 30 31 32]33 34 35 pM [37 38 39 40]41 42 43 44] 46 47 40]49 50 51 52/53 54 55 5657 58 59 60]61 62 63 64]65 66 67 GB]69 70 71 72]73 74 75 76|77 78 79 80 











88 8 8/8 8 8 

1.2.3.4/5 6 7 8/9 10 11 12/13 14 15 16/17 18 19 20/21 22 23 24)25 26 27 28/29 30 31 32/33 34 35 26137 38 39 40/41 42 43 44/45 46 47 45149 50 51 52/53 54 55 56/57 58 59 G8)G1 62 G3 64\65 6G G7 Gs1G9 70 71 72)73 74 75 76|77 78 79 80 

9999/9999199991999 919M 9919995999919 99 919 99.9199991999 51999 9199991899 91899 9199951999 499991999 919998 

12.3415 6 7 B19 1011 12119 14 t5 t6l17 18 19 20lzt 22 23 24125 26 27 26129 30 31 32139 34 35 36137 3B 39 40l4t 42 43 44]45 46 47 4249 50 51 52153 54 50 54197 58 59 GOLGT 62 63 G4l65 G6 G7 Blea 70 71 72173 74 75 76177 78 79 86 
IBM 733727 

i srk a 


5 6 7 89 10 11 12/13 14 15 1} 18 19 20] 18 19 20]21 22 23 24]25 25 fi 26[29 30 31 32]33 34 35 fM]37 36 39 40/41 42 43 44}my 46 47 48]49 50 51 52/53 54 55 56]57 58 59 60]61 62 63 64]65 66 67 GB]69 70 71 712|73 74 75 76|77 18 79 BO 
0 
4 
1 





1888 8/8888 

123415678 
9999/9999/9995/9999/99991999 999999999999 T99GIIDIIIIF OST SSESasisgaggaisggg ggg sy gg sigggggggg 
12.3415 6 7 819 10 11 12113 14 15 16117 18 19 20121 22 23 24,25 26 27 20120 90 31 32139 34 35 36197 3B 39 dolar 42 49 44145 46 47 45149 50 51 SZI52 54 58 SGIST SH $9 GulGT G2 G3 64165 66 G7 Galea 70 71 72173 24 75 26177 78 79 86 


sBM 733727 
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OUTPUT CARDS FOR SLOPE STABILITY ANALYSIS PROGRAM 5171 





fi 2 3 4[5 6 7 8|9 10 11 12]13 14 15 1) 18 19 20]21 22 23 24 25 25 py 28]29 30 31 32[33 34 35 [37 38 39 40/41 42 43 0 46 47 48]49 50 51 52]53 54 55 56/57 58 59 60]6i 62 63 64]65 66 67 68[69 70 71 72]73 74 75 70]77 78 79 80) 








GR EGE CECE EGCG CECEC HEE CECE CCGG EGE EEG 666666666 


12:3.4/5 6 7 8/9 10 41 12/13 14 16 16]17 18 19 29/21 22 23 24)25 26 27 28)29 30 31 32)3y 34 25 36/37 38 39 40141 42 43 44/45 46 47 48/49 59 51 52/53 54 55 56/57 58 59 60)61 62 63 64165 64 67 68189 70 71 72173 74 75 76/77 78 79 20 
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1.2.3 4/5 6 7 8/9 1011 12/13 14 15 16]17 19 19 20121 22 23 24125 26 27 2029 30 31 32193 34 35 90|37 39 39 40141 42 43 44145 45 47 40]19 50 51 52/53 54 55 59/57 5S $9 cvIGT G2 63 G4[65 66 G7 Gsle9 7071 72/73 74 75 ICI77 7B 79 80 
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1.2.3 4]5 6 7 8/9 1011 12/13 14 15 16/17 18 19 20121 22 23 24)25 25 27 20129 30 31 92)33 34 35 36197 39 99 40141 42 49 44145 49 47 42149 50 51 G¥/59 54 55 Sol57 98 59 GUIG1 62 63 64/E5 63 67 EL|G9 70 71 72173 74 75 75}77 78 79 BC 
$999/9999/9999)9999/9 9H 91999 c9999999 919999999 S99 GMs 99 G9 99 O19 e os G9S999II9 G9 899999919998 
1.2.3 415 6 7 819 1011 12113 14 15 16117 19 19 20121 22 23 24,25 25 27 28129 30 31 32133 34 35 36137 39-99 40l41 42 49 44f4s 46 47 48149 50 51 Selsa 54 85 SoIS7 50 G4 GalG1 62 63 G4l65 G6 67 68l69 70 71 72173 74 75 T6l7? 78 79 85 
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SLOPE STABILITY ANALYSIS 
PROGRAM # 5171 


Symbol Table 


Before Instruction No. 212: 


ALY (1) 
AAX(n) 


yee 


=-tAecoord., Of° Lst..end «point, of. .line 
- Y coord. of lst end point of line 
- X coord. of 2nd end point of line 
= Yecoond. of 2nd end point of line 
- Unit weight of material below line 
- Friction coeef. below line 


) - Cohesion coeff. below line 


- X coord. of unordered intersection 
- Y coord. of unordered intersection 
- X coord. of ordered intersection 


- Number of lines in cross-section 

- X coord. of center of failure arc 

- Y coord. of center of failure arc 

- Radius of failure arc 

- Index for counting number of unordered intersections 

- Slope of line 

- (quadratic formula) 

- (quadratic formula) 

- (quadratic formula) 

=“ B - 4(A)(C) | 

- Computed X coord. of simutaneous solution of circle 
and straight line 

- Index for counting number of ordered intersections 

- Y coordinate 

- X coordinate 

- lowest factor of safety for failure arc centers on 
a horizontal line 


- Identification No. 

- Initial radius 

- Degree of consolidation of material below line 

- Index for counting failure arcs having the same 
center 

- Index for counting failure arcs centers 

~ Present lowest factor of safety 

- Index used in ordering intersections 
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After Instruction 


SLOPE STABILITY ANALYSIS 


PROGRAM # 5171 


Symbol Table 


No. 212 (note - variables remaining the same are 


not redefined) 


Summation of cohesion forces 
Summation of tangential or driving forces 
Summation of frictional forces 
Chord length between successive intersections 
on are 
Cohesion force for an increment of arc 
Slice width 
Dist. from origin to centerline of slice 
Summation of normal weight components 
(per arc segment) 
Index for counting number of slices 
(10 per arc segment 
Y distance to bottom of slice on center line 


~ Index for counting number of ordered intersections 


on centerline of slice 


-~ Y coord. of intersection on centerline of slice 
- Y coord. of ordered intersections on centerline 


of slice 

Unit weight of material associated with unordered 
intersections on, centerline 

Y coord. of ordered intersections on centerline 
of slice 

Unit weights associated with ordered intersections 

Total weight of slice 

Summation of frictional forces per segment of arc 

Factor of safety 

Degree of consolidation associated with ordered 
intersections 

Weight of slice for computing normal forces 

Used as index for computed go to 

Unit weight associated with ordered intersections 
on centerline of slice 

Index for counting intersections on centerline 
of slice 

Weight of segment of slice 

Y distance to bottom of slice 

Factor of safety stored for comparison 

Radius of the most critical failure arc at a center 
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SLOPE STABILITY ANALYSIS 


PROGRAM # 5171 


BLOCK DIAGRAM 


OCTOBER 1963 
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SLOPE STABILITY ANALYSIS 


BLOCK DIAGRAM 










READ CARD WITH 





lo. Now ANO Now oF LINES 












READ CARDS WITH COORDINATES 


OF LINES AND PROPERTIES OF SOIL 


READ CARD WITH CENTER COORDINATES 
AND RAOIUS OF STARTING 
FAILURE ARC. 





DETERMINE COORDINATES OF INTERSECTIONS 
OF FAILURE ARC AND LINES MAKING UP 
SECT IONe 


ARRANGE INTERSECTIONS IN ORDER OF 
DECREASING MAGNITUDE OF X=COORDINATE 
AND NUMBER THEM, 









HORIZONTAL DIS= 






ee) ee TANCE BETWEEN TWO SUCCESSIVE YES 

INTERSECTIONS GREATER THAN 
Divide SECTION ABOVE FAILURE ARC AND Divi0E SECTION ABOVE FAILURE ARC AND 
BETWEEN VERTICAL LINES DRAWN THROUGH THE BETWEEN VERTICAL LINES DRAWN THROUGH 
TWO SUCCESSIVE INTERSECTIONS INTO Two (2) THE two (2) SUCCESSIVE INTERSECTIONS 
VERTICAL SLICES OF EQUAL WIDTH. INTO TEN (10) VERTICAL 8LICES OF 


EQUAL WIDTH, 











DETERMINE ANO NUMBER Y=COOROINATES OF 
INTERSECTIONS OF LINES MAKING UP SEC= 
TION WITH THE CENTERLINE OF THE SLICE. 
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SLOPE STABILITY ANALYSIS 
BLOCK DIAGRAM 





To EACH INTERSECTION A88IGN W, Cy, and U OF 
THE SOIL BELOW THE CORRESPONDING LINE. 
To THE SLICE ASSIGN THE F OF THE SOIL 
BELOW THE LOWEST INTERSECTION. 





CoMPUTE THE WEIGHT OF EACH 
SEGMENT OF THE SLICE. 


L 













8&8 THE WEIGHT OF 


YES 
THE SEGMENT NEGATIVE? . 


NO MULTIPLY WEIGHT OF SEGMENT TO BE USED 


1N COMPUTING THE RESISTING FRICTIONAL 
FORCE BY THE DEGREE OF CONSOLIDATION. 


USE ABSOLUTE VALUE OF WEIGHT OF SEG= 
MENT IN COMPUTING THE ORIVING FORCE. 








W/ 





ComMPUTE NORMAL AND ORIVING FORCES 


ON BASE OF SLICE. 





ComMPUTE COHESIVE FORCE ON BASE OF S8LICE BY 
MULTIPLYING COHESION BELOW THE LOWEST 
INTERSECTION ON CENTERLINE OF SLICE BY THE 
CHORD LENGTH OF THE BASE OF THE SLICE, 
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SLOPE STABILITY ANALYSIS 


BLOCK DIAGRAM 


ee OS 









STORE CUMULATIVE TOTALS OF 





NORMAL, ORIVINGsy AND COHESIVE FORCES. 





( 7 \ NO 











4 


18 THIS THE j 
LAST SLICE BETWEEN TWO is 
SUCCESSIVE INTERSECTIONS ON 

THE FAILURE ARC? 
| YES 


FIND FRICTIONAL FORCE ON THIS SEGMENT 






OF FAILURE ARC. 


STORE CUMULATIVE TOTAL OF |: 


FRICTIONAL FORCE. 





LAST SEGMENT OF FAILURE ARC. 


YES 


COMPUTE THE FACTOR OF SAFETY. 
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SLOPE STABILITY ANALYSIS 
BLOCK DIAGRAM 












8 THIS 
THE FIRST FACTOR OF 
SAFETY COMPUTED FOR THIS CENTE 
OF FAILURE ARC? 






YES 







Is THIS FAC= 
TOR OF SAFETY SMALLER 
THAN THE PREVIOUS ONE? 





INCREASE RADIUS BY FIVE FEET. 


PUNCH THE PREVIOUS SAFETY FACTOR, THE 
PREVIOUS FAILURE ARG RADIUS, ANDO THE 
COORDINATES OF THE FAILURE ARC CENTER. 
















HIS FAILURE ARC CENTER 
THE FIRST ONE INVESTIGATED HAV= 
ING THIS Y=COORDINATE? 





A CENTER BEEN IN= 
VESTIGATED WITH A Y=COOROD!INATE 
FIVE FEET LESS? 















ALL CENTERS 
INVESTIGATED HAVING 
‘THIS Y=<COORDINATE, 1S THERE 
ONE AT WHICH THE FACTOR OF SAFETY 18 
LOWER THAN AT CENTERS FIVE FEET 
FROM IT HORIZONTALLY ON 
EACH 81DE7? 











INCREASE X=COORDINATE OF FAILURE ARC 
CENTER BY FIVE FEET. 


RESET RADIUS TO GIVE THE STARTING 
DEPTH OF FAILURE ARC. 





THE LOWEST FACTOR 
OF SAFETY FOUNO OCCUR AT A 
CENTER WITH A Y=COORDINATE FIVE FEET 
LESS THAN THE PRESENT 
Y=COOROINATE? 







YES. 
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SLOPE STABILITY ANALYSIS 
BLOCK DIAGRAM 





SET X=COORDINATE TO GET A FAILURE ARC CENTER 
AT A POINT FIVE FEET FROM THE CENTER WITH THE 
LOWEST FACTOR OF SAFETY ON THIS Y=-COORDINATE 
AND UNOCCUPIED BY ANY OTHER FAII.URE ARC 
CENTER. 


RESET RADIUS TO GIVE THE STARTING DEPTH 


OF FAILURE ARC. 












SET X=COORDINATE EQUAL TO THE 
X-COORODINATE OF THE CENTER WITH 
THE LOWEST FACTOR OF SAFETY. 


RESET RADIUS TO GIVE THE STARTING OEPTH 





OF FAILURE ARC. 
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SLOPE STABILITY ANALYSIS 
BLOCK DOIAGRAM 












THE MINIMUM FACTOR OF YES 


SAFETY AT THIS CENTER LOWER THAN 
AT THE CENTER FIVE FEET 
BELOW IT? 


INCREASE X=COORDINATE BY FIVE FEET. 


RESET RADIUS TO GIVE STARTING 


DEPTH OF FAILURE ARC. 












THE MINIMUM FAC]= 
TOR OF SAFETY AT THIS CEN= 
TER LOWER THAN THAT AT THE CENTER 
WITH BOTH X= ANDO Y=COORDINATES 
FIVE FEET LESS. 
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SLOPE STABILITY ANALYSIS 


PROGRAM # 5171 


FORTRAN LISTING 


OCTOBER 1963 
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C SLOPE STABILITY ANALYSIS WITH SEARCH PATTERN (REVISED) 5171 | 


“DIMENSION X1(10)9X2(10),5Y1(10) 5 Y2(10) »W(10),F (10) »U(10) 





DIMENSION AAX(14),FS(10) 





~ DIMENSTON C (10) pALX014) pATY(14) pATFO1L4) AIC(14) sAIU(14) 





2 READ 45,1Y,IDyNO,N 


Aas SUSAR CEES (B00 GCSES OCC SII IE NPI eSNG RC oe Pe DY 
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SLOPE STABILITY ANALYSIS 
PROGRAM NUMBER 5171 


EXPLANATION OF TERMS 





UNIT GROUPING CARD #1 
Name §—-—s- Field Columns _ Explanation 
Year XX 1-2 Year of Test 
—_— X 3 Minus Punch 
District XX Lh - 5 District Number 
= x 6 Minus Punch 
Number XXX 7-9 Test Sequence Number 
N XX ll 1 Number of Lines in 


Cross Section 
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SLOPE STABILITY ANALYSIS 
PROGRAM . NUMBER 5171 


EXPLANATION OF TERMS 


UNIT GROUPING CARD TYPE 2 
Name Field Columns Explanation 
X1 *XXXX. 1l- 6 LT. Horizontal Coordinate 
: . of Line 
Y1 *XXX. 10 21n LT. Vertical Coordinate 
-of Line 
X2 «XXX. 17 - 21 RT. Horizontal Coordinate 
of Line 
Y2 —*XXX, Pihiceoe RT. Vertical Coordinate 
of Line 
W XXX 30) 93h Weight of Material (P.C.F.) 
F X XXX 37 ~ 41 Friction Angle in Degrees 
C AX. Lh - 48 Cohesion 
U X XXX 51 - 55 Consolidation Ratio 
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SLOPE STABILITY ANALYSIS 
PROGRAM NUMBER 5171 


EXPLANATION OF TERMS 





UNIT GROUPING CARD TYPE 3 
Name Field ~ Columns Explanation 
H XXX, l- 4& Horizontal Coordinate 
of Radius 
G XXX. 6-9 Vertical Coordinate 
of Radius 
R XXX. t] — 1, Radius 
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